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Abstract. We prove large deviation principles (LDP) for the invariant measures of 
the multiclass totally asymmetric simple exclusion process (TASEP) and the multiclass 
Hammersely-Aldous-Diaconis (HAD) process on a torus. The proof is based on a com- 
binatorial representation of the measures in terms of a collapsing procedure introduced 
in [5] for the 2-class TASEP and then generalized in [5], |10) and [11) to the multiclass 
TASEP and the multiclass HAD process. The rate functionals are written in terms of 
variational problems that we solve in the cases of 2-class processes. 



1. Introduction 

In recent years several new and interesting results have been obtained in the study of 
fluctuations of interacting particle systems. Some of these results concern the computa- 
tions of large deviations rate functionals for specific models. 

Given a stochastic interacting particle system, a problem of interest is the determination 
of its invariant measures. When the model is not reversible and the detailed balance does 
not hold, this can be a difficult task. Typical examples are boundary driven stochastic 
lattice gases. Depending on the model you can have available some representations of the 
invariant measures or not. We will give in this introduction a short outline of some of the 
recent progress in computation of the rate functionals of large deviations for the empirical 
measures of the invariant measures in models of this type. The results are interesting for 
several reasons. One reason is that the measures have often long range correlations and 
the corresponding rate functionals are not local. They have a structure very different from 
the one obtained in the case for example of Gibbs measures where you have an integration 
of a function of the density profile. Another reason is that interacting particle systems are 
very effective models of statistical mechanics and the results obtained give insight for the 
behavior of more complex models or real systems. 

From one side there are combinatorial representations for the invariant measures of ex- 
clusion like models starting from which it is possible to compute the corresponding rate 
functionals (see [5] for a recent review and references therein). These combinatorial rep- 
resentations are built up from products of operators satisfying appropriate commutation 
relations. 

A different approach is based on dynamical arguments. Fluctuations of the invariant 
measure can be recovered from fluctuations of paths of the processes. The static rate 
functional for the invariant measure is then obtained solving a variational problem for the 
dynamic rate functional. This leads to an Hamilton-Jacobi equation as a central object (see 
[3] for a recent review and references therein) . Differently from the exact solution approach 
the dynamic one is insensitive to small perturbations of the dynamics corresponding to the 
same macroscopic structure of fluctuations. Nevertheless the Hamilton-Jacobi equation 
is, in general, difficult to solve. 
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In this paper we prove large deviation principles (LDP) for the invariant measures of 
the multiclass totally asymmetric simple exclusion process (TASEP) and the multiclass 
Hammersely-Aldous-Diaconis (HAD) process. Our proof is based on an exact combina- 
torial representation of the invariant measures. Configurations of particles distributed 
according to the invariant measures are constructed applying a deterministic transfor- 
mation, the collapsing procedure, to configurations of particles distributed according to 
product of uniform measures. We generalize the collapsing procedure up to let it act on 
positive measures. An application of the contraction principle allows one to obtain the 
final result. The rate functionals that we obtain are not local. In the case of 2 class models 
they are obtained from a geometric construction on density profiles. 

As our results are obtained from the contraction principle, the rate functionals are 
naturally expressed as infimum of auxiliary functionals. Several non local functionals 
obtained as rate functionals of particle systems are represented in terms of either infimum 
or supremum of auxiliary functionals. An interesting question is whether there is always 
a representation of them as an infimum. More precisely if those rate functionals can be 
obtained from wider LDP using the contraction principles. This is the case of the present 
paper. This is also the case of the invariant measures for the TASEP with boundary 
sources as represented in [7]. The corresponding application of the contraction principle 
has been done in [3]. In [3] it is also suggested that this could be the case for the KMP 
model with boundary sources. A rate functional obtained from contraction of a convex 
rate functional is not necessarily convex, this is also the case of the present result. 

The paper is organized as follows. In section [2] we define the multiclass TASEP and 
the multiclass HAD process constructing them using the basic coupling. In section [3] 
we describe the collapsing procedure as it acts on configurations of particles. We discuss 
also briefly its main properties. In section d] we show how the collapsing procedure is 
used to construct the invariant measures of the processes. In section [5] we generalize the 
collapsing procedure defining its action on positive measures. We discuss also its main 
properties. In section [6] we define empirical measures and relate them to the collapsing 
procedure. In section [7] we derive, from well known results, large deviation principles for 
uniform distributions. In section [8] we derive LDP for the invariant measures of the 2-class 
TASEP and the 2-class HAD process. Using the contraction principle we have the rate 
functionals in a variational form for which we can find the unique minimum in terms of a 
concave envelope construction. These are theorems 18.21 and 18.51 that are the main results 
of the paper. The rate functionals are not convex. We discuss also in variational terms 
the typical density of first class particles when the total density is known and the typical 
total density when the density of first class particles in known. In section [9] we derive 
LDP for the multiclass TASEP and the multiclass HAD process. The rate functionals 
are written in terms of variational problems. It is interesting to study such problems and 
others proposed in the section. We obtain also a recursive relation. 

The combinatorial constructions of [2], [9], [TO], [llj are different from the original 
solution of the 2-class TASEP obtained in [6], that is based on products of non commuting 
operators. It is interesting to derive the same result of the present paper starting from this 
alternative solution. A probably different representation of the same rate functional will 
maybe appear. It is also interesting to study the problem using the dynamic approach. 
The major problem here is the lack of a complete dynamical LDP. A problem of interest 
is also the study of the variational problem (j9.ip for models with more than 2 classes of 
particles. 
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To avoid confusion we remark that wc use similar symbols for mathematical objects 
that play a similar role in the TASEP, the HAD process and in the general framework 
of positive measures. We use also the same symbol C to indicate the collapsing operator 
both when it acts on configurations of particles and on positive measures. This is due to 
the fact that the second one is a natural generalization of the previous one. 



2.1. TASEP. The totally asymmetric simple exclusion process (TASEP) is a model of 
stochastic jumping particles satisfying an exclusion rule. Let Ztv := Z/{NZ) be the 
discrete one dimensional torus with N sites. Every site x € Zjv can be either empty or 
occupied by a particle. The state space of the process is Xj^f := {0, 1}^^. Given t] G 
a configuration of particles, we will say that the site x G Z^r is occupied by a particle if 
r]{x) = 1 and empty otherwise. Every particle waits an exponential time of rate one and 
then tries to jump to its nearest neighbor site to the left. If the site is already occupied by 
another particle then the jump is suppressed. This informal description can be summarized 
by the following generator of the dynamics 



where Cx is the oriented bond {x, x + 1) and rj^"^ is the configuration of particles obtained 
from r] rearranging the values at the extremes of ex in decreasing order, according to the 
following definition 



A generalization of the previous model is obtained labeling some of the particles as first 
class particles and the remaining ones as second class particles. When a first class particle 
tries to jump over a second class particle it succeeds and the two particles exchange their 
positions. When a second class particle tries to jump over a first class particle the jump 
is suppressed. The natural state space for such a process is {0,1,2}^^ obtained from 
the choice of assigning value to empty sites, value 1 to sites occupied by a first class 
particle and value 2 to sites occupied by second class particles. We will instead describe a 
configuration with a pair (771,772) with both 771 and 772 elements of Xj^. The configuration 
7/1 G X]\f is such that 771 (x) = 1 when in x there is a first class particle and 771 (x) = 
otherwise. The configuration 772 G X^ is such that 772 (x) = 1 if in x there is either a first 
or a second class particle and 772 (x) = otherwise. We endow Xn with the natural partial 
order -< defined from 



2. MULTICLASS MODELS 



^iv/(77) := Yl U{v'^)-f{v)) 




77^^ ^7/(x)^^(x), yxeZN- 

By definition rji ^ 772 so that (771, 772) G where 

■= {{Vi, ■■■,Vk) -Vi^ Xn, rji < m+i} ■ 



The two descriptions of the state space are equivalent and a bijection between 
{0, 1, 2}^^ is defined from 



and 



^(x):=inf{z:77i(x) = l, z = 0, 1, 2} , 
where we defined 770 (x) := 1 — 772 (x). 
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The generator of the above described 2-class TASEP is 

LNf{vi,m)-= {fir]T,r]2"') - f{vi,V2)) ■ 

This generator clearly defines also a jointly, order preserving, evolution of two TASEP, 
usually called basic coupling. 

A further natural generalization, called the k-class TASEP is obtained introducing par- 
ticles of class up to a fixed natural number k ^ N. When a particle of class i tries to jump 
over a particles of class j with j > i then the positions of the two particles are exchanged. 
When a particle of class i tries to jump over a particle of class j with j ^ i the jump 
is suppressed. The state space is now I^^ that is in bijection with {0, 1, ... , /c}^^. The 
generator of the dynamics is 

Ln f ivi, Vk) ■■= Yl {fivl",---,rik')- f{Vi,---,Vk)) ■ 

This generator clearly defines also a jointly, order preserving, evolution of k TASEP, 
usually called basic coupling. 

2.2. HAD process. Let A := M/Z be the one dimensional torus. The Hammersley- 
Aldous-Diaconis (HAD) process (after [T2] and [1]) is a stochastic evolution on finite 
subsets of A. Let 

Qn ■= {x:= {xi, . . . ,xn} ■■ XieA;Xi^Xj i j} 

be the collection of all finite subsets of A with points and let ^1 = Dn^n. Labels to 
points are given in such a way that Xj+i is the nearest point of x to the right of x,. Given 
X € an initial condition, the HAD process preserves the number of points and is defined 
as follows. Every point Xj waits an exponential time of rate |(xj,Xj+i]| and then jumps to 
a point uniformly chosen in (xj, Xj+i]. This dynamics can be easily summarized from the 
following generator 

Lf{x) = f du (/(x«) - /(x)) , 

where the set x" is defined from 

u j Xi if u ^ (xj,Xj+i] , 

'~ \ u if ue {xi,Xi+i] . 

This formula holds for any i = 1, . . . , |x|, with the convention X|a,|_|_i := xi. With proba- 
bility one all the points stay distinct along the evolution. 

The multiclass HAD process has not a simple and intuitive behavior. As in the case of 
the TASEP a natural way to define it is through the basic coupling. The state space is 
^k,-,N, defined as 

x»GJ7^,; x«Cx(^+i)} . 

Here ^ A'^i ^ A''2 • • • ^ are k natural numbers. This is a natural set to describe 
points in A with associated an integer class from 1 to k. Points are the elements of x^'^). 

The class associated to x^*^^ is inf |j : x^*^^ G x^-'-'l. The multiclass dynamics is defined 

from the following generator 

Lf{x^^\...,x^^^) = j du (/(xW'",...,x('^)'") -/(xW,...,x('=))) 
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that describes also a joint, inclusion preserving, evolution of k HAD processes. This joint 
evolution is usually called basic coupling. 

Both for the TASEP and the HAD process a natural way to introduce basic coupling is 
via a graphical construction where all the coupled processes evolve using the same random 
marks. In the case of the HAD process for example these marks are points of a rate one 
Poisson point process on the cylinder A x See [9], [TO], [11] for a more detailed 

description. 



To describe the invariant measures of the multiclass processes previously introduced we 
have to explain a collapsing procedure introduced in [2j. We start describing its action on 
configurations of the TASEP. 



The collapsed configuration [^i] is obtained from a mass preserving (i.e. number of 
particle preserving) transformation of the configuration tji . The transformation is defined 
algorithmically as follows. Give any order to the particles of the configuration r]i. Move 
the first particle of this configuration that is on a site, say x, such that r]2{x) = 0, to 
the first site on the right, say y, such that: r]2{y) = 1 and r]i{y) = 0. Update the rji 
configuration according to this movement and iterate the procedure using at every step 
the same order fixed at the beginning. The final configuration does not depend on the 
specific order chosen. 

We briefly discuss some of the properties of this algorithmic transformation. For more 
details we refer the reader to |2j where this construction was introduced, and to [9], [TO], 
|llj for an interpretation in terms of queue theory. 

By definition we have that Cr,2[77i] ^ r/2. Let E{y,x) be the excess of r/i particles in 
[y, x] defined as follows 



ze[y,x] 

Lemma 3.1. There is a positive flux of particles across the bond if and only if there 
exists y E Z^r such that E{y,x) > 0. 

Proof. If y is such that E{y, x) > then in [y, x] there are more particles of iji type than 
of r]2- During the collapsing procedure some particles will necessarily flow out of [y, x] and 
this can happens only through Cx- Conversely let us suppose that for any y € Z^r we have 
E{y, x) ^ 0. Let us order the particles of iji from right to left starting from x and let 
Zi be the site corresponding to the i particle. Due to the fact that E{zi,x) ^ 0, particle 
number 1 will be allocated on a site belonging to [a^i, x]. Due to the fact that E{z2,x) ^ 0, 
particle number 2 will be allocated on a site belonging to [^2,2;] (different from the one of 
particle number 1) and so on. No particles will flow across Cx- □ 



3. Collapsing particles 



Let us call 



Given ^ Mi ^ M2 ^ two natural numbers, we define a collapsing operator 

that maps the pair of configurations (r/i, 7^2) into the pair of configurations 

C(r/i,r?2) := (C^2 [r/i], r?2) . 





(3.1) 
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Let [•]_,_ be the positive part defined as 



X if X , 
otherwise 



Lemma 3.2. The total flux of particles across is 

J{x) ■= sup [E{y,x)]j^ . 
y&N 

Proof. When J{x) = this follows directly from the previous lemma. If J{x) > we can 
argue as follows. Clearly J{x) is a lower bound of the total flux, because for any y the 
excess of rji particles E{y,x), if positive, has necessarily to flow across ex- Let y* be the 
first element of Z^r to the left of x such that E{y*,x) = J{x) > 0. There is no flux of 
particles across ey*-i. This follows from the previous lemma and the fact that for any z 
it holds E{z, y* - 1) ^ 0. In fact, if z G [x + 1, y* - 1] then 

E{z,x) ^ J{x) = E{y\x) , 

that implies E{z, y* — 1) = E{z, x) — E{y*,x) ^ 0. If instead z E [y*, we have 

E{z, y* + E{y\x) = E{z, x) + E{x + 1, x) 

^ J{x) + Ml - M2 = Eiy*,x) + M1-M2 , 

that implies E{z, y* — 1) ^ Mi — M2 ^ 0. We established that all the particles flowing 
across ex were originally in [y* , x]. From the characterization of y* we deduce immediately 
that E{y*,z) > for any z G Let us order the particles of 772 type contained in 

[y*,a;] from left to right and let Zi be the site corresponding to the particle number i. 
Remember that there is an excess J{x) of rji particles in this interval. Due to the fact 
that E{y* , zi) > a particle of type tji will be allocated in zi. Due to the fact that 
E{y* , Z2) > a particle of r/i type will be allocated in Z2 and so on. At the end all the 
sites Zi will be occupied by rji particles and exactly the excess of particles J{x) will flow 
trough Cx- □ 

Lemma 3.3. For any interval [a, b] it holds 

E Cr,Am]{^) = E ^i(^) + -^(^ - 1) - '^(^) • (3.2) 

x£[a,b] x£[a,b] 



Proof. This property follows directly from the conservation of mass. Equation (13. 2p simply 
states that the number of r/i type particles that are at the end of the collapsing procedure 
in the interval [a, b] is obtained from the number of particles present initially plus the 
number of particles entered from the left side minus the number of particles exit from the 
right side. □ 

The collapsing procedure is defined in a similar way for configurations of the HAD 
process. Let ^ A^i ^ A''2 be two integer numbers. We define the collapsing operator 

T 

that maps the pair of configurations (x, y) into the pair of configurations 

C(x,y) := (^Cy[x],yj . 

The collapsed configuration Cy[x} is obtained moving to the right points of x. The trans- 
formation is defined algorithmically as follows. Give any order to points of x. Move the 
first point of x that does not belong to y to the nearest point of y to the right that does 
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not belong to x. Update the x configuration according to the previous transformation 
and iterate the procedure. The final configuration does not depend on the specific order 
chosen. 

Also in this case it is possible to define an excess E{u, v) oi x points in the interval [u, v] 
E{u,v) := {xr\[u,v]} — {yn[n, u]} 
and consequently a flux of x particles at v 

J{v) := sup[E{u,v)]_^_ . 
ueA 

Note that J is right continuous i.e. it holds J{v) = J{v~^) := lim^iQ J{v + e). All the 
lemmas previously listed, appropriately reformulated, holds also in this case. We do not go 
into details here because the collapsing procedure will be generalized to a wider framework 
in section [5l We only write down the analogous of equation ()3.2p in this case 

jCj^y n [n,t>]| = {xn[n,u]} + J{u~) - J{v) , (3.3) 

where J{u^) is the left limit of J at n, i.e. J{u~) := lim^iQ J{u — e). 

4. Invariant measures 

4.1. TASEP. For the TASEP the number of particles Ylxez^^ ri{x) is a conserved quantity. 
For any fixed integer M ^ N, the TASEP with M particles is an irreducible finite state 
Markov chain on Xi^^m and consequently has a unique invariant measure. The process has 
then a one parameter (i.e. M) family of invariant measures that is easily seen to coincide 
with the family of uniform distributions of M particles on Z^r 

j,fj{n) = \ (m) ^ if ri(^ Xn,m , 
I otherwise . 

These are all the extremal invariant measures; all the remaining are obtained as convex 
combinations. 

The 2-class TASEP conserves the number of first class particles X^ojgZjv 
number of second class particles '^x^i^{'r]2{x) — rii{x)). For any fixed pair of non negative 
integer numbers Ai and A2 such that Mi := Ai and M2 := Ai + A2 ^ N, the TASEP 
with Ai first class particles and A2 second class particles is an irreducible finite state 
Markov chain and has an unique invariant measure. 

The result in [2j states that this invariant measure is 



<ix<^joC-^ (4.1) 

We use the symbol x to indicate the product of measures. In general given a measure 
and a measurable map T with the symbol /i o T^^ we denote the pull-back measure 
defined from 

(/.oT-i) (A) :=f,{T-HA)) 

for any measurable set A. 

We remark that (j4.ip is a measure on iff . This two parameters (Mi and M2) family 
of invariant measures constitutes all the extremal invariant measures. 

To give a combinatorial representation of the invariant measures of the /c-class TASEP 
we need to extend the collapsing procedure of section [3l This extension to the case of 
more than two classes of particles is contained in [9] and further discussed in [10], pT] . 
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Let Ai, . . . , Afc be k non negative integer numbers such that X^^Lj^ Aj ^ N. Let us call 
also Mj := J2i=i We define a collapsing operator 

Cfc : Xn,Mi X Xn^M2 X • • • X X^^Mk -^AT^ 

that associates to the configurations (771, . . . , 77^,) the configurations 

defined as follows. The configuration S,k coincides with r/^. The configuration £,k-i coincides 
with Crj,.[r]k-i]- The configuration coincides with C^^. [C^j._ Jr/fc_2]]. In general the 
configuration S,k-j is obtained from the composition of j collapsing procedures 

e^-, := C,, [C,,_, [. . . C,,_^^, . . . ]] . (4.2) 

Obviously according to this definition C = C2. The result contained in [9j states that the 
invariant measure for the TASEP with Aj i-class particles (i = 1, . . . , /c) is 



This k-parameter family of invariant measures constitutes all the extremal invariant mea- 
sures. 

4.2. HAD. The HAD process conserves the number of points. In $7 tv the unique invariant 
measure is given by the support of the values of i.i.d random variables uniform in A. 
Equivalently the invariant measure is a uniform Poisson point process in A, conditioned 
to have N points. These are all the extremal invariant measures, all the remaining are 
obtained as convex combinations. 

The 2-class HAD process conserves the number of first class points and the number of 
second class points. Let Ai and A2 be two non negative integer numbers and let Ni := Ai 
and A2 := Ai + A2. The 2-class HAD process with Ai points of first class and A2 points 
of second class has a unique invariant measure that has a combinatorial representation 
in terms of the collapsing operator. In fact the result in [11] states that this invariant 
measure is 

These are all the extremal invariant measures, all the remaining are obtained as convex 
combinations. 

The generalization of this combinatorial construction to the case of more than two 
classes of points is described also in fll]. We proceed as in the case of the TASEP. Let 
Ai,A2,...Afc be non negative integer numbers and call Nj := X]i=i^«- We define a 
collapsing operator 

Cfc : 0,Ni X X ■ ■ ■ X ilTVfe ^Ni,...,Nk 

that associates to the configurations [x^^\ ■ ■ ■ ,a;('^)) the configurations 

defined as follows. The configuration y^^^ coincides with x^''\ The configuration y^*-'"^) 
coincides with C^(fc) [x^*''^^)]. The configuration yC^"^) coincides with C^(k) [C^(k-i) [x^'^"^)]]. 

In general the configuration yC^"-?) is obtained from the composition of j collapsing 
procedures 



(fc-i) 



...a 



(fc-j+i) 



(4.3) 



FROM COMBINATORICS TO LARGE DEVIATIONS 



9 



The result in [TT] states that the invariant measure for the HAD process with Aj i-class 
particles (i = 1, . . . , A;) is 

These are all the extremal invariant measures, all the remaining are obtained as convex 
combinations. 

5. Collapsing measures 

We start with some definitions. The set of positive measures on A will be denoted as 
^A. With M.^ we will denote the subset of Ai of measures p absolutely continuous with 
respect to Lebesgue measure and with M^'^ the subset of Ai^ containing the elements 
such that their densities satisfy the condition 

dp 

^ ^ 1 a.e. . 

du 

With abuse of notation we will indicate with p both a generic element of A4^ and its 
density. Finally, given m G M+, we call Mm := {p & A4 : Jj^ dp = m}. Likewise we set 

Ml := A^o n Mrn and M^' := f] Mm- 
We define a partial order ^ on A4 saying that 

Pi^ P2 <^ dpii^ dp2 
J A J A 

for any measurable ^ C A. We then call 

/'''^ := {{pi, . . . ,pk) ■■ Pi e M; Pi ^ pi+i} . 

We want to generalize to this framework the algorithmic constructions illustrated in 
section [3l We are dealing no more with configurations of particles, but with positive 
measures. Given pi G Mmi and p2 G Mm2 with mi ^ m2 we want to define a collapsing 
operator 

C : Mm, X Mm2 ^ I^'^ 

that associates to the pair {pi,p2) the pair 

C(/9l,P2) := {Cp2[pl\,P2) , 

where the collapsed measure Cp^ [pi] ■< p2 is obtained moving mass of pi to the right. The 
natural way to define such a procedure is through a generalization of equations ()3.2p and 
()3.3p . We start defining the excess of mass E{u, v) of the measure pi in the interval [u, v\ 
as 



E{u,v) := 




By definition E{u, v) is right continuous in v (i.e. E{u, v^) := lim^io E{u,v + e) = E{u, v)) 
and left continuous in u (i.e. E(u~,v) = E{u,v)) . It satisfies also some simple addition 
rules 

E{a,c) =E{a,h-) + E{h,c) 
E{a,c) = E{a,b) +E{b+,c) 

In the above formulas we use the convention E{u, u~) 
the fiux of mass across v € A defined as 

.J{v) := sup [E{u, v)] 



V 6 G [a, c] , 

V 6 G [a, c] . 

= E{u^ , u) := 0. Then we introduce 
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We give the following definition of the collapsing operator. 

Definition 5.1. The collapsed measure Cpjipi] is defined in such a way that for any 
interval (a, b] it holds 



dCpM = / dpi + J(a) - J(6) . (5.1) 

{a,b] J(a,b\ 

Note that the action on intervals of this type completely defines the measure. 
Let 

J ■.= {veK: 3n G A s.t. E{u,v) > 0} . (5.2) 

Given v (z J^, let u such that E{u, v) > 0. From the right continuity in v we have 
E{u,v^) = E{u,v) > and consequently there exists an e > such that v + 6 £ J' for 
any ^ 6 ^ e. From this we can deduce that J' = UiJi. Where Ji are at most countable 
many disjoint intervals either of the type \li,ri) or {li,ri). This fact can be proved with 
an argument very similar to the one used to characterize open sets on M (see for example 
section 11.11 of [E]). The condition mi ^ m2 implies the fact that the strict inclusion 
i7 C A holds. A sketch of the proof is as follows. Let, by contradiction, assume that 
J' = A. Then for any w G A you can prove there exist an interval whose right closed 
boundary is v, it is either open or closed at the left boundary and is such that 



< J{v) = dpi- dp2 ■ 

Moreover any other interval with the same property is contained inside U^. Then for any 
vi V2 we have that either Uy-^ R = or they are one contained inside the other. We 
define Wv ■ ^{weA-. v&Jw}^^- Given vi ^ V2 then either = '^^ they are disjoint. 
Moreover 



dpi - I dp2> G A 

SO that at most countable different Wv can exist and they form a partition of A. We finally 
have 

^ mi - m2 = X\ dpi - Jj^ dp2 

= Ev, (/w., dpi - J^^^ dp2) > , 

a contradiction. An alternative route is obtained showing with arguments similar to the 
ones of the next lemma that there exists an element of A in a neighborhood of the left 
boundary of any that does not belong to J'. 

Lemma 5.2. Consider v G J^i. If JTi = [li,ri) then 

J{v) = E{k,v) . 

If instead Ji = {li,ri) then 

J{v)=E{lt,v) . 

Proof. Here and hereafter in some proofs we need to distinguish the two cases when J7i = 
[li,ri) or J'i = {li,ri). We will give the proofs in the case Ji, = [li,ri). The proofs for the 
other case are analogous. Given v € J we consider Wn a maximizing sequence for J{v), 
i.e. a sequence such that 

lim E{wn,v) = J{v) . 
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If Wn G {v,li) then 

E{wn,v) = E{wn, 1-) + E{k,v) ^ E{k,v) , (5.3) 
where we used the definition of and the fact that there exists an element of J^'^, the 
complement of in A, in any neighborhood of li. Without loss of generality we can then 
consider Wn G [h , v] . Note also that the simple argument in (j5.3p also implies that for any 
V ^ Ji there exists u € \li-,v\ such that E{u,v) > 0. 

We prove now that for any w € J^i it holds E{li^w) ^ 0. Let us introduce 

N ■.= {w(^Ji: E{k,w) < 0} . 

Due to the fact that E(u, v) is right continuous in v we have that Af is the union of at most 
countable many intervals either of the type [aj,bj) or of the type {aj,bj). We give the 
proof in the case there are only a finite number of intervals and oi, the nearest boundary 
element to the right of k, is the boundary of an interval of the type (ai,6i). As before, 
the proof for the remaining cases is analogous. As already showed there exists u £ [li,ai] 
with E{u, ai) > 0. Then we have 

^ E{k,af) = E{li, ai) = E{li,u-) + E{u, ai) ^ E{u, oi) , 

a contradiction. This imply that TV = 0. Finally we get for the maximizing sequence 
Wn £ [k , v] that 

E{li,v) = E{li,w:^) +E{wn,v) ^ E{wn,v) . 
This means that Wn := k is a maximizing sequence and this implies the first statement of 
the lemma. □ 

Note that a direct consequence of this lemma is that J is right continuous, an important 
fact to have that definition 15.11 is well posed. 

It is possible to introduce a measure 7 defined from 



dj := J{b) - J{a) . 

'{a,b] 

The measure 7 is not a positive measure. In fact we have 



d-f = lim / dj = lim {J{v) - J{v + e)) = . (5.4) 



40 J{v+t,v] 

Definition (15. ip then becomes 

Cp2 [Pl] = Pl-1 

and using (j5.4p we derive the conservation of mass 



/ dCp^lpi] = dpi- ^7 = / dpi ■ 
Ja J a J a J a 



We give now a simple representation of the collapsed measure Cp^lpi]. Similar represen- 
tations hold also for collapsed configurations of particles in the TASEP and points in the 
HAD process but their generalization to the case of positive measures is not straightfor- 
ward. 

Lemma 5.3. The collapsed measure Cp^[pi\ has the following representation 

Cp2 [pl] = PiXjc + p2Xj + X] '^P^ ~ '^P^ ^"-^ ' (^•^) 

where with the symbol Xa denote the characteristic function of the set A and with 5^ 
the delta measure in v. 
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Proof. We show that the weight given to any interval (a, b] from the measure defined by 
the right hand side of (15. ip coincide with the weight given to the same interval from the 
measure on the right hand side of ()5.5p . This implies that the two measures coincide. 
We need to verify the following identity 



/ dpi + J(a) - J(6) 

J{a,b] 

I dpi+ I dp2+ ( I dpi - I 

J{a,b]nj'= J{a,b]nj r,t(a,b] ^ 



This is equivalent to 



dpi - I dp2 = - \ dpi - dp2] Xj (a) 

{a,b]nj J{a,b]nj y[h*M J [I,*, a] J 

+ (/ dpi- ( dp2\xj{b)+ il '^P^- I '^P^) ■ (^-6) 



In the above formula, when a G ^ we call Ji* the interval to which it belongs and we 
assume it is of the type [/i*,rj»), when 6 G ^7 we call Jj* the interval to which it belongs 
and we assume it is of the type [lj*,rj*). The proof in the remaining cases is similar. Note 
that it is possible to have i* = j* . 

The sum on the right hand side of (j5.6p can be written in the following way. 



riG{a,b] ^-^^^ -^^^ ^ 

U dpi- I dp^ + U dpi- I dp^xAa). (5.7) 



Using this formula the right hand side of (15. 6p becomes 



JiC(a,6] 



V ( / fipl - / c?P2 ) + / dpi- dp2 \ Xj{a) 

+ I / dpi- dp2] Xjib) =y2{ '^P^ ~ / '^P'^ 

'[lj,,b] J[lj*M J i \JJ,n{a,b] JJrn(a,b] 



that clearly coincides with the left hand side of (15. 6p . □ 

We have the following simple properties. 
Lemma 5.4. The measure Cp^ [pi] is a positive measure. 

Proof. We use the representation (j5.5p . The first two terms on the right hand side are 
clearly positive. In the case Ji = we have that 



/ dpi - f dp2 = E{h,r-)^0 , 



as shown during the proof of lemma [5^21 This implies that also the third term is positive. 
An analogous argument holds in the case Ji = {li,ri). □ 
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Lemma 5.5. It holds 

CpM<p2. (5.8) 
Proof. Using formula ()5.5|) we obtain that condition ()5.8p is equivalent to 

PiXjc - P2Xjc + ^ {^j^ dpi - dp2^ Sr.^O. (5.9) 

We show that the weight associated to any interval (a, b] from the measure on the left 
hand side of ()5.9p is negative and this implies the statement of the lemma. 
We need to prove that for any (a, 6] it holds 



/ dpi- dp2+ I dpi- dp2] ^0 



(5.10) 



We use the fact that 



dpi - / dp2 
dpi- - / dpi - / dp2 

{a,b] J{a,b] ^ \JJ^n{a,b] J J,n{a,b] ^ 

/ dpi- dp2- [ dpi- dp2 ] 

J{a,b] J{a,b] J,C{a,b] "^^^ ' 



dpi - I dp2 \ Xj{a) - / dpi 



f dp2]xj{b) . (5.11) 



In the above formulas we use for the labels i* and j* and the corresponding intervals, the 
same convention as in lemma [5^31 Using (15. IIP and equation ()5.7p we obtain that the left 
hand side of ()5.10p is equal to 



dpi - dp2+ \ dpi- dp2 Xj {a) 

(a,b\ J{a,b] \J[h*,a] J[k,,a] J 

dpi- [ dp2]xj{b). (5.12) 

[i^,,b] J[i^,,b\ ) 

In the case i* 7^ j* we have that (j5.12p can be written as 

E{k* , lJ*)Xj {a)Xj {b) + E{li* , h)xj {a)Xjc (b) 

+E{a+ , lj.)xjc {a)xj (b) + E{a+, b)Xj. {a)Xjc {b) 

and all the terms are non positive. 

The case i* = j* can happen only when a £ J' and h £ J . In this case we have that 
(j5.12p is equal either to or to 



dpi - / dp2 ^ . 

□ 

When both pi and p2 belong to M.^ then E{u, v) is continuous in u and v and conse- 
quently E{li,r~) = E{li,ri). From what we proved in lemma [5?2] we have that E{li,r~) ^ 0. 
From the fact that G J''^ we have that E{li,ri) ^ 0. This implies that E{li,ri) = 0. 
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Then the third term on the right hand side of (15. 5p is not present and in this case we 
simply have 

Cp2[Pl] = PlXjc +P2Xj ■ 

This means that also Cp^ [pi] € with a density a.e. given by 

CMiv) = \ l % ^^1' (5.13) 



P2 4 



The collapsing operator C is not continuous with respect to the weak topology. This 

(n) 

can be easily seen from the following example. Consider the sequences of measures p\ := 

(n) 

^fL-Ll\ Py :=5i+53. Clearly we have 

V2"*"nJ 2 4 

Moreover it holds 

(n) n^oo , (n) n— >oo 

p\' ^ pi and P2 ^ P2 
with pi := 5i, p2 ■= Si + Ss and convergence is in the weak topology. We have also 

2 2 4 

Cp2[pl] = ■ 

As a consequence 

hm C(pi"\p(")) / C(pi,p2) = C ( hm 

n— »oo \n— *oo 

Nevertheless it holds the following continuity result. 

Lemma 5.6. The collapsing operator C is continuous with respect to the weak topology 
on M'i^ xTW^^. 

Proof. Let {p^i^\p'^'^) G x M.m2 be a sequence of measures weakly converging to 

(/5i)P2) G -^mi X ■^%i2- Then for any fixed u G A we define the following nondecreasing 
continuous functions on the interval [0, 1] 

![v-uM '^P^ ^f""^ ^) ' 

mi if u = 1 , 



G2(u) :-- 



I[v-u,v] dp2 ifue [0, 1) , 
m2 if u = 1 . 



Likewise using the measures p^^ and pj"^ , we define the functions G^"^ and Gg"^ that are 
nondecreasing but not necessarily continuous. The weak convergence implies the pointwise 
convergence of G^"^ to Gi and G'^^ to G2. Monotonicity of all the G, and continuity of 

(n) 

the limit functions imply also the uniform convergence. Consider for example G\ and 
Gi. Fix an arbitrary e, let 



Ui := inf {u : Gi(n) = ze} ; i = 0, . . . 
and define also ^i|^iinj_)_i := 1- For any i let n{i) be such that for any n > it holds 
\G^\ui) — Gi{ui)\ ^ e. Let also n* := maxj {n(i)} . Monotonicity says that when u G 

[Ui,Ui+l) 

G^r\ur) ^ g'CHu) ^ G^;'\u,+,) 
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and 

Gi{ui) ^ Gi{u) ^ Gi{ui+i) . 

These inequalities imply 

\G^r\u) - Gi{u)\ ^ m^x[\G^^\u,)-Gi{ui+i)\,\G^;'\u,+i)-Gi{ui)\] . 

When n > n* both terms inside the max are ^ 2e. Let us show this for example for the 
second one. It holds 

\G^^\u^+l) - Gliu^)\ ^ \Gt\u,+i) - Gi{ui+i)\ + \Gi{u,+i) - Gi{ui)\ 

and both terms on the right hand side are ^ e. This implies the uniform convergence of 
G^"^ to Gi. The uniform convergence of G^"^ and G'^^^ implies the uniform convergence 
of 



G^^^ — Gg"^ to [Gi — G2]_^_- This implies the convergence of 



j(«)(i;)= sup 

ue[o,i] '- 



G^r\u)-Gr{n) 



to 



J{v) = sup [Gi{u) - G2{u)]_ 

ne[0,l] 



(5.14) 



(5.15) 



for any v (z A. To prove this simple fact we see that, calling u* a maximum point in (j5.15p . 
liminf„^oo J^^-^v) = liminf„^oo ( sup^gro ii G^j^\u) - 03"^ (n) 



> liminf^ 



G'i'\u*) -G)^''{u*) =J{v). 



-40, 



On the opposite direction, let us call w'^ a maximizing sequence in (|5.14p . i.e. a sequence 
such that 

Gi")«)-G(")«)" . 



lim 

fc— ►oo 



j(")(v 

Consider now a sequence k{n) such that 

jW(t;)-fGi")(u;^(„))-G(")«(„)) 



1 

< - . 

n 



It holds 



hm sup J(") (v) = hm sup G^^^ (t.^(„) ) - G^^^ (t.^(„) ) 

n^oo n— >oo L 

We then obtain 

limsup„^oo J^^\v) = lim sup. 

= limsup„_oo 



1 y'^Hn)! ~ ^2 y'^kin)' 

Gi«(„))-G2«(„))j^J(.). 



We conclude the proof showing that 



(n) 
Pi 



lim„_oo rfpi"^ + J(")(a) - 

I{a,b] dpi + J{a)- J{b) = J^a,b] dCp, [pi] 
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The convergence of f^^ dp^^ to JJ.^ dpi follows from the weak convergence of pS"^ to pi. 

The pointwise convergence of J(") to J has been shown above. This implies the statement 
of the lemma. □ 

Also in this framework the collapsing operators Ck are defined likewise in ()4.2p and 

()4.3p in the cases of particles configurations. The continuity properties of are derived 
also from lemma ISTHl 

6. Empirical measures 

Given a configuration rj Xj^ we associate to it its empirical measure iTNiv)- This is 
an element of Ai defined as 



N 

Given a collection (rji, . . . ,rjk) of /c configurations of Xf^f we will write 

7r^(ryi, . . . ,ryfc) := (vrAr(?7i), . . . ,TTN{Vk)) ■ 

It is easy to check that 

■^N iC^2[Vl]) = C7rjv(r,2)[^Jv(^l)] • (6-1) 

This is an important identity that justifies the fact that we used the same symbols for the 
collapsing operators in different frameworks. Note in fact that in (|6.1|) we are using the 
same symbol C with different meanings. Equation (|6.1|) is the key identity to check the 
following commutation property 

TT^ o C = C o vr^ . (6.2) 
The validity of (j6.2p directly implies its generalization 

7T%oCk = CkOTT% . (6.3) 

Note that the same commutation relations hold also in the case we had given a slightly 
different definition of empirical measure associated to a configuration of particles of the 
TASEP. Sometimes the empirical measure associated to a configuration r] € X]\[ is defined 
as an element of A4^'^ whose density is a.e. 



Also in this case dSl]), ([621) and I^Ml hold. 

In the case of the HAD process there is not a natural scale parameter as in the case 
of the TASEP where it is the size of the lattice. We will consider families, with index 
a natural number N, of HAD models containing [Nm] particles, with m a positive real 
number and [•] the integer part. The scale parameter is and for any configuration xgQ 
we define the empirical measure as 

i 

Note that we do not require in this definition that coincides with \x\. As in the case of 
the TASEP, given a collection {x^^\ . . . jX^'^^) of k configurations of we will write 

Also in this framework the analogous of (16. ip . (16. 2p and (16. Sp hold. 
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7. Large deviations for uniform distributions 

In this section we quickly derive, from well known results, large deviations principles 
for the invariant measures of the TASEP and the HAD process with particles of only one 
class. 

Given a sequence pN of probability measures on a Polish metric space X we say that it 
satisfies a large deviation principle (LDP) with parameter and rate function / : X — > 
M+ U {+00}, if 

limsupA^"Mog/U7v(C) ^ - inf I{x) VC C A , C closed ; 

V^+00 a^GC 

limmf N-^logfiN{0) ^ - inf I{x) VO C A , O open . 

V— >+oo x£0 

The rate function / is lower semicontinuous and is called good if it has compact level sets. 

Let us recall a classical LDP result for sampling without replacement as stated for 
example in [4j. 

Let y^^) = |y[^\ • • • , I be A^ elements of A such that 

N 



1=1 

where — > is the weak convergence and A € . Let be the uniform measure on M-uples 
of elements of y^^\ M ^ A^: 



1/^ 



'•••'^«M /; \ otherwise . 



This is the measure obtained from a procedure of M samplings without replacement among 
the A^ elements of y^^^ ■ Consider now the case M = AI(N) = [A^m], where [ • ] is the 

integer part and m € (0, 1]. Then when A^ diverges the measures v^^^^ o vr^-*^ satisfy a 
LDP on A4 equipped with the weak topology with parameter A^ and with a good and 
convex rate function given by 



mH 



a) +(i-m)F(^|A) zf X-peMi-m, (7^^) 



-00 otherwise 

Where -ff (•!•) is the relative entropy. 

To get the rate functional for the invariant measure of the TASEP with [Nm] particles 

we need to consider y^^^ = j^. With this choice the measure A coincides with the Lebesgue 
measure and consequently 

/a Piu) log p{u)du if peM" , 
+00 otherwise . 



H{p\X) 



In this case the condition for finiteness in ()7.ip is equivalent to /3 G Alm^. Observing that 
in this case the collection of all the M-uples of y^^^ is in bijection with Aat^jv/ (using the 
correspondence rj(x) = 1 if and only if belongs to the M-uple), we obtain the following 
result. 
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Proposition 7.1. Let vfj he the invariant measure of the TASEP with M particles. When 

N diverges the family of measures z/jy^™' o -k'^ satisfy a LDP on A4 equipped with weak 
topology, with parameter N and with good and convex rate function 

Slip) •= I -^^^^ hm{p{u)) if M^yi , 
\ +00 otherwise , 

l-x 



where hm{x) = x log ^ + (1 - x) log j 

We use the symbol 5*1 to remark that it is the rate function for a one class system. The 
parameter m is understood. 

A LDP for the invariant measures of the HAD process with particles of only one class 
follows directly from Sanov theorem (see [3] section 6.2). Sanov theorem states that if Xi 
are i.i.d random variables taking values in A and having common law A then the empirical 
measure 

1 ^ 

i 

satisfies a LDP on M equipped with the weak topology with parameter and with a 
good convex rate functional given by the relative entropy 

H{p\X) if peMi, 
+00 otherwise . 

From this general fact we can easily deduce the following proposition as a special case 
when A is the Lebesgue measure. 

Proposition 7.2. Let pM be the invariant measure of the HAD process with M points 
and let m be a positive real number. When N diverges the family of measures P[Nm] satisfy 
a large deviation principle on Ad equipped with the weak topology with parameter N and 
with the good and convex rate functional 

S (p) ■= / -^A^^ kmipiu)) if pe , 
^ I +QO otherwise , 

where km{x) = xlog^. 

For simplicity we use the same symbol Si already used for the TASEP, but the rate 
functions are different. 

8. LDP FOR 2-CLASS PROCESSES 

Theorem 17.11 is immediately generalized to the case of product measures. Consider 
< mi ^ 7712 ^ 1 real numbers. The family of measures 

.ir"^^lx.!,^-]lofvri,)-^ (8.1) 



satisfy a LDP on A^ x 7W endowed with weak topology, with parameter A^ and with good 
and convex rate function given by 

TA=ii'Kduhm,{pi{u)) if p^eM'^A^^, .g2) 
+00 otherwise . 

We are interested in proving a LDP for the empirical measures of the invariant measures of 
the 2-class TASEP. This means that we are interested in proving a LDP for the sequence 
of measures 

x.!,^-l)o(C2)-^]o(.^)-^ , (8.3) 
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that due to identity (16. 2p coincides with the sequence of measures 

o(C2)"^ . (8.4) 

Lemma 15.61 suggests that we can apply the contraction principle. 

The contraction principle (see for example [3] section 4.2.1) states that if /U at is a se- 
quence of measures satisfying a large deviation principle on a Polish metric space X with 
a good rate functional I{x) and f : X \s a, continuous map from X to another Polish 
metric space Y , then also the sequence of measures /x^r o satisfy a LDP on Y with 
good rate functional K given by K{y) = infj^jg^-i^y)}. /(x). This formulation can in fact 
be extended (see remark c at page 127 of [4J) to the case when the map / is continuous 
only at the elements x G X such that I{x) < +oo. This is exactly our setting. 

Before stating the theorem that we obtain following this strategy we need some facts 
and notations. Consider an interval U := [u', u'"] C A and a density p{u). Let us introduce 
the extended function 

pu^^ ._ f Jul p{w)dw if ueU , 

P ^ ' ' y — oo otherwise . 

Note that in the case of measures in Ad^ we can use the above notation for integration 
because there is no difference between open and closed intervals. Let us also consider the 
extended function on the real line M defined as 




F]^{u^ + u) if ue [0,\[u\u'-]\] , 
— oo otherwise . 



We call the concave envelope of F^(n). We then define the following extended function 
on A 

Fjfiu) :=¥^^{\W,u]\) 

Note that the following properties hold 

Fjf in) ;,Fjf(u) , 
Fjf{n^) = Fjf{u^) = 0, 
F^{un = F^{un = i:: p{w)dw. 

Finally we call p^ a density on A such that p^{w)dw = Fjf{u) for any u ^ U and 
p^ {u)xi^c{u) = a.e.. It can be shown that p^ is a positive measure and it belongs to 
when p G M^'^. 

Given pi and p2 measurable functions on A, the seiU := {u G A : pi{u) = P2{u)} is a.e. 
equivalent to the disjoint union of at most countable many closed intervals Hi := \u\, u[] . 
Indeed it is simple to check that the family of sets with this property is a cr-algebra that 
includes all the open sets. Hence it must include all the Borel sets. 

Lemma 8.1. Let (pi,p2) € /^'^ n (A^^i ^ -^m2)- ^'^^'^ ^^'^^ ^^'^ P'^^'^ iipi,i^2) £ 

^ -^m2 ^-^ such that C{ipi,Tp2) = {pi, P2) if and only if the following conditions are 
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satisfied 

Proof. We recall that due to the fact that the measures are absolutely continuous with 
respect to Lebesgue measure, the above functions F are continuous on the interior part 
of the intervals where they are different from — oo. First we show that given a pair 
(V'i,V'2) £ -^mi ^ -^m2 ^^^^ satisfies conditions (|8.5p then we have C(^i,'i/'2) = {pi,P2)- 
Clearly we need only to prove that Cp2[^i] = pi. Remember that for measures belonging 
to the collapsing procedure acts as in (I5.13p . Conditions number three and four in 
(18. 5p imply that for any u (zUi and for any i we have 



dipi - d/)2 ^ . 

J u 

Let us consider the following subset of U 

i 

From the fact that equality between the F holds at the boundary of each Ui and from 
continuity in the interior part we deduce that this is an open set. We show now that U 
coincides with the set J' as defined in equation (15. 2p for the pair of measures (^1,^2)- 
Clearly U C. J . This follows from the fact that for any u gU with u €Ui it holds 



/ dV'i - / 

Jul Jul 



dp2 > 



To prove that U = J we need to show that for any u (zU^ and for any f € A it holds 

pu 

dtpi - dp2^0 . 



Note that 

dipi - dp2 ^0 , 

[v,u]nu'' J[v,u]nU'^' 

due to the fact that on we have a.e. that ipi{w) = pi{w) and pi{w) < P2{w). We 
consider first the case u gW^. In this case we have 



#1 - I dp2 = I dipi- dp2 



+ 5Z f/^^i-/ dp^ + ( f ^' dil,i- f dp^Xuiv) ■ (8.6) 

When V (z U we called Ui* the interval to which it belongs. All terms on the right hand 
side of (|8.6p are nonpositive. Consider now the case u G lADlA^ and call lAj* the interval to 
which it belongs. In this case we need to modify formula (j8.6p multiplying the last term 
on the right hand side by X[„ {u^j*) and adding 



d^i - I dp2] X[„,„]c(nj.) + 
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Both terms are nonpositive. Note that i* can coincide with j*. 

Conversely we assume that for (V'ii'02) £ -^m^ x ■^m2 holds C{'ipi,'ip2) = {pi,P2) and 
show that this implies the validity of conditions (j8.5p . 

The validity of the first condition in (jS.Sp is obvious due to the fact that the collapsing 
operator C preserves the second component. 

Using (I5.13P and the previous statement we have that a.e. it holds the following equality 

Cpiii^iKu) = 'il^i{u)xjc{u) +p2{u)xj{u) = pi{u) . 

We multiply both sides by x^c (u) and obtain 

(Mu) - Pi{u))Xjcnuciu) + ip2{u) - piiu))Xjr,uc{u) =0 , a.e.. (8.7) 

Due to the fact that the two terms can be different from zero on disjoint sets we have that 
is equivalent to the two equations 



(ipiiu) - pi{u))xjcnu'^iu) = , a.e. , ( 

{p2iu) - pi{u))Xjnuciu) =0 , a.e.. ( 
From the fact that pi{u) < P2{u) a.e. on and from equation (j8.9p we deduce that 



d\= dX , (8.10) 

where A is the Lebesgue measure. Equation (jS.Sp imposes that tpi{u) = pi{u) a.e. on 
J''^ n and using (j8.10p we obtain the validity of the second condition in (j8.5p . 

In the case of absolutely continuous measures we have that J is an open set and this 
implies that u\ and u'^ do not belong to J" for any i. In fact let us suppose for example 
that m[ G J', then there exists a ball Bu^{e) centered in such that Bu^{e) C J' and 
Sb r(e)nw<:'^'^ ^ ^ (recall that A denotes Lebesgue measure). This is impossible due to 

i 

(|8.10p . From the fact that u\ and u[ belong to J''^ we deduce that 

J{uD = J{ul) = . 

Using this we obtain 

dp2 = / dCp,[il^i] = dVi , (8.11) 



that is the third condition in (j8.5p . 

Let us suppose that there exists an u € Z^j such that F^^{u) < F^^{u). Then using also 
()8.1ip we obtain 



dil^i- / dp2 = F^;{u)-F'll{u)>Q . 



This implies u\ ^ J that is impossible. As a consequence we obtain the validity of the 
fourth condition in (18.51). □ 



Now we can state and prove our large deviations result. 

Theorem 8.2. Let < mi < m2 ^ 1 real numbers. Consider the 2-class TASEP on 
Ztv with respectively \Nmi\ first class particles and [Nm2] total particles. When the 
pair {r]i,r]2) is distributed according to the invariant measure of the process, we have that 
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7rjj[r]i,r]2) satisfies a LDP with parameter N and with good rate function S2{pi, P2) de- 
fined as follows. It takes the value +00 when {pi,p2) l'^'^ H ^A^m^ x A^m2^- When 

{pi,P2) £ l'^'^ n ^A^mi X A1m2^ it takes the value 

S2{pi,P2) = / hrm{pi{u))du + 'y] / hmAp^'iu)) + / hjn2{p2{u))du , (8.12) 

where Ui are disjoint closed intervals andU = UiUi is a.e. equivalent to the set 

{li G A : pii^i) = P2{u)} . 
The symbol p^' is defined before the statement of lemma \8.1\ 

Proof. As outlined before we can apply the contraction principle. The sequence of mea- 
sures in (jS.ip satisfies a LDP with good rate functional given by ()8.2p and, as shown in 
lemma 15.61 the map C is continuous on every point where ()8.2p is different from +00. 
We immediately get that the sequence of measures (|8.4p satisfies a LDP with a good rate 
given by 



S2{pi,P2)= inf ( / hm^{il;i{u))du + / hm2iip2iu))du 

(8.13) 

By convention the infimum over an empty set is defined as +00 and remember that we 
are calling ipi both the measures and the corresponding densities. From ()8.13p we see 

immediately that S2 is equal to +00 on (^I^'^ n ^A^m^ x A1m2^^ • 

The case mi = m2 = m is not covered by the theorem (in fact it could be, choosing 
appropriately an interval lA). In this case the rate functional is different from +00 only 
when pi = p2 and consequently U = A. The rate functional is then 

S2(pl,P2) = I L^rn{p{u))du if p^ = p2 = p (£ M'^A^ , 

\ +00 otherwise . 

When mi < m2 we have a strict inclusion U C A. 

First we prove existence of a minimizer for the variational problem (I8.13P , then we prove 
uniqueness and finally we characterize it. 

In the case of the TASEP, existence of a minimizer can be derived directly showing that 
()8.13p is a minimization problem of a lower semicontinuous functional over a compact set. 

We show instead a more involved proof that works also in the case of the HAD process. 
When {pi,p2) G l"^'^ then C{pi,p2) = {pi,P2)- This means that the we can modify the 
infimum in ()8.13p restricting to pairs (^1,-02) that satisfy the further condition 



(V'i,V'2) G 



2 2 

-M^' X M^'^ : ^ / hmMi^))du ^ E / hmMN)du \ . (8.14) 
i=i i=i ''^ ) 



Due to the fact that (18. 2p is a good rate function we have that ()8.14p is a compact set. The 
constraints set in ()8.13p is easily seen to be a closed subset of TWmi x (and this holds 

also without the condition of bounded densities) . In fact consider a sequence (V'l"'* , V'2"'* ) 
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belonging to this set and converging to {i^i,tp2) ^ A^mi x A^mi- Then by lemma \5M we 
have that 

C(Vi, V2) = C f lim = hm €(^^1"^ V?^) = {pu P2) , 

\n—>oo / n^oo 

which means that the set is closed. We obtained an infimum of a lower semicontinuous 
functional over a compact set and the existence of a minimizer follows. 

We prove now uniqueness of the minimizer. The functional to be minimized in (18.130 
is strictly convex. This follows directly from the fact that the real functions h^yi^ are 
strictly convex. The set on which we are minimizing is also a convex subset of A^mi x 
Aim2- This follows directly from its characterization given in lemma [8Tl If {'ipi,ip2) and 
(V'i,V'2) satisfy conditions (18. 5p then clearly also the convex combination (^ic,'02c) := 
c{tjji,'ip2) + (1 — , -02) satisfies the same conditions and this is clearly true also for the 
additional condition of bounded density. A classical result in convex analysis [8] guarantees 
uniqueness of the minimizer of a strictly convex functional over a convex set. 

Finally we characterize the unique minimizer. 

Using conditions ()8.5p we write ()8.13p as 

S2{Pl,P2) = Jf,hm2{P2{u))du + hmiiPliu))du 

where 

:= G M'^' : F^^ ^ F^- F^^{u^) = Fjf^{ul)] . (8.15) 
We need then to study the variational problems 



inf / hmi{(f)i{u))du . (8.16) 

Existence and uniqueness of the minimizers can be shown as before. It remains to char- 
acterize them. 

From the strict convexity in x of the real function hm{x) and Jensen inequality we have 
for any interval [^1,1)2] Q A and any density (j) 

77 77 / hm{(piu))du ^ hm 77 77 / (j){u)du . (8.17) 

Moreover this inequality is strict as soon as the density (piu) is not a.e. constant. A 
geometric interpretation of this inequality is the following. Let (p and (p' be two densities 
defined on an interval V C A and let [fi,i'2] ^ V. Consider the case in which F^{ u) = 
F^,{u) for any u {vi,V2) and the graph of F^, when u G [fi,f2] linearly interpolates 
{vi,F^{vi)) and (^2, ^^(^2)). More precisely 

F^'iu) = ^^^'''y^f.^'''\ \[v^,u]\)+F^{v^) ; u G [v^,V2] . (8.18) 

Note that if G A4^'^ then necessarily also (p' G A4^'^. Inequality (j8.17p then simply says 
that 

'' hrrMu))du ^ [ hm{<p'{u))du , (8.19) 



with the strict inequality holding if (p{u) and (p'{u) do not coincide a.e.. In the rest of 
the proof we consider pairs {<p,(p') whose corresponding F^^{u) and F'^^{u) are related as 
before. We can then apply inequality (|8.19p . 
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Consider a S Ai, such that F^' is not a concave function. Then clearly there exists 

a (/)' such that F^,''{u) ^ F^^{u) for any n. In particular this implies that also 0' G Aj. 
Inequality (j8.19p implies that (j) can not be the unique minimizer. The unique minimizer 
(f)* of (j8.16p has then necessarily F^l concave. 

Let us now consider (j),ip ^ Ai such that both F^' and F^* are concave and moreover 
there exists an u £h(i such that F^'(n) > F^'(u). Then can not be the unique minimizer. 
We can in fact construct a (p' considering an affine function through (ti,F^'(ti)), obtained 
using an element in the superdifferential of F^' at u. We then have 



F^; (v) ^ min { F^^ {v),F^^{v)} , V^; G Ui 



so that (j)' G Ai. Inequality ()8.19p then implies that <j) can not be the unique minimizer. 
The unique minimizer <p* is then necessarily such that is the smallest among all the 

concave functions that are above Fp», that is its concave envelope F^^. This shows that 
on the interval Ui we have 0*(n) = (}^^{u) = p^^{u) a.e.. □ 

The rate functional 5*2 is non negative and is zero if and only if ipi{u) = nii a.e. and 
tp2iu) = m2 a.e. and corresponding p2iu) = m2 a.e. and pi{u) = Cm2[w-i]('u) = w-i a.e.. 

The rate functional 5*2 is not convex. This can be shown by the following example. Let 
us consider densities {pi,p2) and [p\, P2) defined a.e. as 

Pi{u) = Xn iXu) ; p\{u) = \x.i Au) ] P2{u) = pI{u) = x.i Au) . 
We consider the convex combination 

{plc,P2c) ■■= c{pl,P2) + (1 - c){pI,P*2) . 

We have that 

S2{pi,P2)= 1 h3.{p2{u))du+ I hi[p^?''^\u)\du+ I hi{pi{u))du, 
[0 — ] 

where p^ (u) = ^x^^ We have also 



S2{p*i,P2)= I h-A{p2{u))du+ / hi{pl{u))du 

and for any c G (0, 1) 

S2{pic,P2c)= I hi{p2{u))du+ I hi{pic{u))du. 



Convexity of S2 would imply the validity for any c G [0, 1] of the following inequality 

cS2{pi,P2) + {l-c)S2{pl,p*2)-S2{pic,P2c)^0 . (8.20) 

If we take the limit c | 1, on the left hand side of (j8.2Up we obtain 

hi{p^i'^\u))du — / hi{pi{u))du , 

that is clearly strictly negative. This shows that S2 is not convex. 
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To simplify notations in the following lemmas the fact that all the measures are abso- 
lutely continuous, have bounded densities and have a fixed total mass will be understood. 
It is understood also the fact that {pi,p2) G -f^'^- 

We can still obtain interesting results from the contraction principle. The following 
identity has to be satisfied 



inf S'2(yOl,P2) = Si{p2) = / hm2{p2{u))du . 

Ja 



121) 



From the microscopic point of view this identity simply derives from the fact that if 
we forget the labels first and second class and we just look at positions of particles the 
dynamics that we observe is a TASEP. From the variational point of view we have the 
following lemma. 

Lemma 8.3. The unique minimizer /)* in (j8.2ip such that 

S2{p*i,p2) = Si{p2) 

is given by 

pl = Cp,[7m], (8.22) 
where mi denotes the measure with a constant density equal to mi. 

Proof. We have 

inf{pi} ™f{(V,i,^2): C{Vl,t^2)=(pi,P2)} (/a {i^l{u))du + hrn2 {ij2{u))du) 

= jf^hm2{p2{u))du + 'm.i^^ fj^hm^{ipi{u))du 
= j\hm2{p2{u))du . 

The unique minimizer in the expression above has a density ipi{u) = mi a.e. and corre- 
spondingly the unique minimizer p^ in (j8.2ip is obtained from the collapsing procedure 
applied to ^J, that is ([02]). □ 

The minimizer p^ in ()8.22p corresponds to the typical density of first class particles 
when the system is conditioned to have a total density p2. 

Still from the contraction principle we have that the following identity has to be satisfied 

inf 5*2(^1,^2) = 5'i(pi) = / hm^{pi{u))du . (8.23) 

From the microscopic point of view this identity simply derives from the fact that if we 
forget second class particles and observes only first class particles what we see is a TASEP. 
Identity (j8.23p can be deduced also from purely variational arguments. 

Fix the density pi and call V = UjVi a subset of A a.e. equivalent to the subset 
{it G A : pi{u) > 7712}. The sets Vj are disjoint closed intervals, Vj := [i^', To any such 
an interval we associate an element w\ K. This element is the nearest if G A to the left 
of v\ such that 

[ \m2 - pi{z)) dz = . (8.24) 

J w 

It exists due to the fact that (m2 — pi{z)) dz = m2 — mi > 0, the function of w given 
Iw ("^2 — Pi{z)) dz is continuous and condition (|8.24p identify a closed set. Consider 
now the intervals Vi := 
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Lemma 8.4. The unique minimizer P2 in ()8.23p such that 

S2{pl,p*2) = Si{pi) , 

has a density a.e. equal to 

p*,{u) = l ^/^^U.V., (8 25) 

[ Pi{u) tf ue UiVi . 

Proof. We have 

S2{Pl,P2) - Si{pi) = f^hm2{P2{u))du + hm^{pi{u))du 
+ [p'i'iu)) - {pi{u)) 

= !u- hm2ip2iu))du + Ju, (pliu)) 
+ Ei lu, {pf{u)) - Ei ^mi (Pl(^)) • 

Where we used the fact that pi{u)Xi,{u) = p2{u)x^{u) a.e.. We add and subtract 
Ei lu ^"^2 ( Z'^' ) ^"^ ^'^'^ finally we obtain 



S2{PI, P2) - Si{pi) = f h^,{p2{u))du + Y^ f hm^ipfiu))^^. (8.26) 

We used the fact that given an interval V and a density p 

{hrriiipiu)) - hmiipiu))) du = ( [ p{u)du] log ^^^^^^ + log 



m2(l — mi) 1 — 777-2 

depends only on mi, m2 and the total mass Jy p{u)du. This imply that 

(hmiiPliu)) - hmi{pi{u)))du= / (^hm2iPi"iu)) - hmj_{p^"iu)fj du , 



due to the fact that pi{u)du = p^'^{u)du. The right hand side of (j8.26p can be zero 



if and only if p2{u)Xuc{u) = m2Xuc{u) a.e. and f}^' {u)Xu.{'^) — '^2Xu X'^) ^^^^ ^"^Y ^• 
This happens if and only if p2 is constructed as in (j8.25p . Let us show this fact. 

Given two intervals Vj and Vj then they are either disjoint or one contained inside the 
other. This follows from the following statement: if w\ then w\ G (f[,u;j). To 

prove the statement observe that, by definition of w\, for any u G {w\,vl) we have 



{1712 - pi{z)) dz + j {m2 - pi{z)) dz = 



and moreover the second integral is strictly negative. As a consequence for any u € {w\,vl) 
it holds 

ru 

/ {m2 - Pliz)) dz > . (8.27) 

If wl [Vj,vl) then J \ {1112 — Pi{z)) dz > 0. Clearly we have also J^^ {m2 — pi{z)) dz < 
for any u G {Wj,Vj). As a consequence we deduce that there exists an n G {w[,Vj) such 
that J^J (?7T,2 — Pi{z)) dz = and this implies the above statement. 

We consider the subfamily of intervals Vi^, composed by the intervals Vj that are not 
contained inside intervals Vj with j i. Note that UfcVjj. = UiVi. 
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The fact that (18.250 is a minimizer of ()8.26p and consequently also of (18.230 follows from 

the fact that inequality (|8.27|) is equivalent to inequality ^ Fpi'' with the equality 

sign holding only at the boundary of Vij, . This implies the fact that F^^'' coincides with 

and consequently p^''^ iu)x~ (u) = ^2X~ (u) a.e.. The fact that (|8.25|) is the unique 

minimizer can be shown from the fact that (18.260 is strictly positive for different density 
profiles. □ 

The density profile P2 in (|8.25p is the typical total density profile when the system is 
conditioned to have a density profile pi of first class particles. 

Following the steps of all the proofs presented for the TASEP you see that the only 
properties of the model that we used are: the strict convexity in x of hm{x) and the the 
fact that hmix) ^ with equality if and only if x = m. Both properties hold also for 
km{x). As a consequence all the above statements hold also for the the HAD process. 
Starting from theorem 17.21 and proceeding as before we obtain the following result. 

Theorem 8.5. Let < mi < m2 positive real numbers. Consider the 2-class HAD 
process on A having respectively [Nmi\ first class particles and [Nm2\ total particles. 
When the pair {x^'^\a^'^^) is distributed according to the invariant measure of the process, 
we have that t^%{3^^\x^'^'^) satisfies a large deviation principle with parameter N and with 
good rate function S2{pi,P2) defined as follows. It takes the value +oo when {pi,p2) ^ 
/^'^ n X A^^a)- When {pi,P2) G -^^'^ n {M'^^ x TW^J ^"^""^ 

S2{pi,P2) = / kmi{pi{u))du + y2 / kmi(pi"{u))+ / km2{p2{u))du , (8.28) 
Jw^ i JUi ^ ^ Jk 

where lAi are disjoint closed intervals andlA = UiUi is a.e. equivalent to the set 

{u £ A: pi{u) = P2{u)} . 

The symbol p^' is defined before the statement of lemma \8.1[ 

The rate functional ^2 is non negative and is zero if and only if ?/'i(u) = mi a.e. and 
■02 (li) = m2 a.e. and corresponding p2{u) = m2 a.e. and pi{u) = (7^2 = mi o-e-- 
Also in this case the rate functional 5*2 is not convex. 

The typical density p^ of first class particles when you condition the system to have a 
total density p2 is given by (|8.22|) . 

The typical total density profile for the system conditioned to have a density profile 
pi of first class particles is given by (j8.25|) . 

9. LDP FOR MULTICLASS PROCESSES 

Theorem 17.11 is also immediately generalized to the case of product of k measures. 
Consider < mi ^ m2 ^ ^ . . . ^ m^ ^ 1 real numbers. The family of measures 

.r^lx...x.r^l)o(vr^^-' 

satisfies a LDP with parameter N and with good rate functional given by 

Yli=l Ia '^mAPii'^)) if Pi e > 

+00 otherwise . 
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We are interested in proving a LDP for the empirical measures of the invariant measures of 
the k-class TASEP. This means that we are interested in proving a LDP for the sequence 
of measures 



[Nmi 



X U 



[Nmk 
N 



(Cfc 



-1 



-1 



that due to identity (|6.2p coincides with the sequence of measures 



[Nrrii 



X V 



TV 



As in the previous section we can apply the contraction principle obtaining the following 
result. 

Theorem 9.1. Let < mi ^ . . . ^ ^1 he real numbers. Consider the k-class TASEP 
on Zat with [Nrrii] particles of class ^ i. When (t/i,...,^) is distributed according to 
the invariant measure of the process, we have that 7r^(7/i, . . . , 77^) satisfies a LDP with 
parameter N and good rate function Sk{pi, • • • , Pk) given by 



Sk{pi 



,Pk) 



inf 




S^{lpl,...,1pk)- ^i&M^^\, Ckilpl,--;1pk) = {pl,-,Pk)^ 

with the convention that the infimum over an empty set is +00. 



hmMiu))du] , (9.1) 



Remember that in (19. ip we are indicating with both the measure and the corre- 
sponding density. The functional Sk is nonnegative and zero if and only if pi{u) = a.e., 



moreover it takes the value +00 on 



'mi 



• • • X M^J', 



We obtained an interesting variational problem that will not be studied in this paper. 
We remark only the following fact. Existence of a minimizer for ()9.ip can be proved using 
the same strategy as in theorem 18.21 Uniqueness of the minimizer is not guaranteed as in 
the case of 2-class models. We have in fact that when k > 2 the set 



{(V"!, • • • ,^fc) : ^i&M^], Cfc(V'i,...V'fc) = (/Oi,---,/0fc)} 



(9.2) 



is not necessarily a convex set. This follows from the following example. Let e be any 
positive real number < | and consider the measures defined from the following densities 



and 



■02 

-01 
V'2 



+ 8xri 1 



4Xri 1- 

U ' 4 



^2 '2 



2xr5 5_^^, , 

4Xr3 3 +4xr3 3 

4X[i l+el + 8X[i 1 



We have that 



C3(01,V'2,03) = C3(0i,V'2,03) 

where the measures pi have densities a.e. equal to 

Pi 



{Pl,P2,P3) , 



P2 
P3 



4X[i 1 

U'4 

4X[i 1 

U'4 

4X[i,i 



+ ^^[|'l + f] 
+ ^^[|'H^1 
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It is easy to check that for the convex combination 

{<t>l,<t>2,<p3) = ^(V'1,V'2,'03) + ^(V'1,'02,V'3) 

we have €3(^1, 02, 03 ) / (pi,P2,P3)- 

The same kind of result is easily derived also for the HAD process starting from theorem 
EH 

Theorem 9.2. Let < mi ^ ... ^ mi be real numbers. Consider the l-class HAD 
process on A with [Nmi] points of class ^ i. When {x^^\ . . . is distributed according 

to the invariant measure of the process, we have that vr^(x^^-*, . . . , x^'^) satisfies a LDP 
with parameter N and with good rate function Si{pi, . . . , pi) given by 

I 



Si{pi,. . . ,pi) = inf , / kmMiiu))du] (9.3) 

with the convention that the infimum over an empty set is +00. 

The functional Si is nonnegative and zero if and only if pi{u) = mi a.e.. Moreover it 
takes the value +00 on [/''^^ n (A^^^ x • • • x M'^^)^. 

Still using the contraction principle we obtain the following identity, valid both for the 
TASEP and the HAD process, whose study from the variational point of view seems to 
be interesting 

inf Skipi, ...,pk)= Sk-i{pi, ...,pj,...,pk) . (9.4) 

{pj} 

With the symbol pj we indicate the fact that the measure pj is missing. From the mi- 
croscopic point of view identity ()9.4p derives from the fact that if you change the class 
of j-class particles to j + 1 the dynamics that you observe is the one of a (fc — l)-class 
process. 

We derive now a recursive relation. We derive it for the TASEP but it holds also for 
the HAD process. To simplify notations the fact that all the measures pi, tpi and (pi 
involved are absolutely continuous, have bounded densities and have fixed total mass will 
be understood. Also the fact that (pi, . . . , pk) G l''"'^ is understood. We can write (19. ip as 

Sk{pi,...,Pk) = fj^hm,,{pk{u))du 
+ inf{(v.i,...,v.,-i): c,(^^,...,^,.„p,)={pi,-,Pk)} (Ei=i lAhrnMiu))du^ , (9.5) 

because if Ckiipi, ■ ■ ■ ,ipk) = {pi,---,Pk) then Vfc = Pk- Let us call {(pi, . . . , (pk-i) := 
Ck-i{ipi, . . . , ipk-i)- We have then that (0i, . . . , (pk-i) ^ ^pk^^ where 

:= {(0i,...,0fc_i) G/^-i'^ C,,[(f>i]=pi, i = l,...,k-l] . 

Then equation (19. 5p becomes 

Skipi, ■■■ ,Pk) = ff^ hmk{pk{u))du 

+ ^"f{(01,...,?ifc_l)GA^-l}i'^f{(^lv..,V'fc-l): Cfe_i(V'l,...,^fc-i) = (9il,...,</'fc-i)} (Etl /A^™.(V'i(^^))^^'") 

that finally becomes the following recursive relation 



Skipi, ■ ■ ■ ,Pk) = / h.mk{Pk{u))du + inf S'fc_i(0i, . . . , 0fc-i) • 

{{0i,...,0fc_i)eA^-i} 
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